arXiv:1503.00160vl [hep-ex] 28 Feb 2015 


PROCEEDINGS 

0F SCIENCE 



Search for CP violation in D° 7i~7i + 7i° decays at 
LHCb 


Shanzhen Chen* 1 ^ 

University of Manchester 

E-mail: shanzhen.chen@cern.ch 

The LHCb experiment has recorded the world’s largest sample of charmed meson decays. This 
paper presents a study of a D° meson decaying into a final state containing a neutral pion in 
LHCb. The search for CP violation exploits a novel model-independent unbinned technique to 
assign a p- value for the no CP violation hypothesis. With a data sample size exceeding that of 
previous measurements by almost an order of magnitude the world’s best sensitivity is obtained. 
The /7-value of no CP violation hypothesis given data of D° —)> k~k + tl° decay analysed is found 
to be (2.6 ±0.5) x lO” 2 . 


Flavorful Ways to New Physics - FWNP, 

28-31 October 2014 

Freudenstadt - Lauterbad, Germany 


* Speaker. 

^On behalf of the LHCb collaboration 


© Copyright owned by the author(s) under the terms of the Creative Co mm ons Attribution-NonCommercial-ShareAlike Licence. 


http://pos.sissa.it/ 








Search for CP violation in D° —)> n 7t+7t° decays at LHCh 


Shanzhen Chen 


1. Introduction 

The D° 7i~7i + 7i° decay 1 is a singly Cabibbo-suppressed neutral D meson decay. The amplitudes 
of this decay channel proceed through c —>> udd , which lead to the final states p°7T°, p + 7T _ , and 
p _ 7T + . These final states are common for both D° and D°. 

Time-integrated CP asymmetries in D° decays can have three components: direct CP violation 
in decays to specific states, indirect CP violation in D° — D° mixing, and indirect CP violation in 
interference of decays with and without mixing. Within the Standard Model, CP violation expected 
in D decay is of O(10 -3 ) or less, however contributions from particles that are not described in the 
Standard Model (SM) and participate in the loops of the penguin amplitude can enhance the CP 
violation effects expected within the SM [1]. Therefore CP violation in D° —»• 7T _ 7r + 7r° decay 
provides sensitivity to physics beyond the Standard Model. 

The previous most sensitive studies of the D° —>> 7r _ 7r + 7r° decay has been done by the BaBar [2] 
and Belle [3] collaborations. At LHCb, with the world’s largest dataset, a novel method, known as 
energy test [4, 5], is used for the first time to search for CP violation in this decay, and the world’s 
best sensitivity of this decay mode is achieved [6]. 


2. Data 

Neutral pions reconstructed in LHCb can be grouped into two categories according to the photon 
cluster structures recorded in the electromagnetic calorimeter (. ECAL ). There are pions that have 
lower momentum (typically p^ < 2 GeV/c) for which both final state photons are reconstructed 
separately in the ECAL (resolved pions), as well as pions that have higher momentum and thus a 
merged cluster reconstructed in the ECAL (merged pions). The two 7T° samples provide coverage of 
complementary regions of the D° —> 7i~n + n° Dalitz phase space and thus the use of both categories 
of pions contributes significantly to the sensitivity of the analysis. 

The data sample used in this analysis is recorded during the 2012 run at a centre-of-mass en¬ 
ergy of 8 TeV, with the integrated luminosity L = 2fb _1 . The offline selections of the two categories 
are performed separately. The resolution of merged pions is worse as there is only one cluster in 
ECAL , nevertheless due to large transverse energy and low combinatorial backgrounds, a better 
purity in merged sample is achieved. After selection, 416 x 10 3 resolved candidates and 247 x 10 3 
merged candidates are obtained with purity of 82% and 91%, respectively. Dalitz plots of resolved 
and merged samples, and the combination of the two are shown in Fig. 1. 


3. Energy Test 


The energy test is an unbinned model-independent statistical method to search for time-integrated 
CP violation in multibody decays. This method relies on the comparison of the probability density 
functions of the two flavor conjugated samples in phase space. A test statistic T is formed to realize 
the comparison, defined as 


i n i n i n,n 
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Figure 1: Dalitz plot of the (a) resolved, (b) merged and (c) combined D° —» 7t 7 r + ^° data sample. En¬ 
hanced event densities in the phase-space corners originate from the p(770) resonances. 


The metric function i///y = yt(dij) = e~ d ^^ 2(J ^ is chosen to be Gaussian function with a tunable 
metric parameter <7, and dij = | Ax z j | = | (m^ 7 — , m^ 7 — , m 2l J — m 2 a ) | is the distance between 

two events in phase space, where the a, b, c subscripts indicate the final-state particles. All three 
invariant masses are used in the calculation of distances between events. Using all three invariant 
masses does not add information, but it avoids an arbitrary choice that could impact the sensitivity 
of the method to different CP violation scenarios. 

The three terms in the T formula sum over the weighted distances among D° decay events, 
among D° decay events, and between events in different samples, respectively. For a CP violating 
sample, events would be more likely to gather at different regions in the phase space for CP conju¬ 
gated decays. Therefore, the average distance in each flavour tagged sample would be smaller than 
the average distance between the two different sets. Since the metric iff (dij) is generally falling 
with increasing distance, the CP asymmetry would lead to a large T value. 

A permutation method is also used to simulate samples without CP violation by randomly 
reassigning the flavour of each candidate. By comparing the nominal T value observed in data to a 
distribution of T values obtained from permutation samples, a p-value under the hypothesis of CP 
symmetry is obtained as the fraction of permutation T values greater than the nominal T value. 

The statistical uncertainty of the p- value is obtained as a binomial standard deviation, reflect¬ 
ing the limit number of permutation tests. If large CP violation is observed, the nominal T value 
tends to be lying outside the bulk of permutation T values (see Fig. 2(a)). In this case, the permuta¬ 
tion T distribution can be fitted with a generalised extreme value (GEV) function, which is defined 
as 



(3.2) 


with normalisation N, location parameter /i, scale parameter <5, and shape parameter £,. The p- value 
from the fitted T distribution can be calculated as the fraction of the integral of the function above 
the nominal T value. This approach is used for the sensitivity studies, which will be described in 
Sec. 4. The uncertainty on the fitted p-v alue is obtained by randomly resampling the fit parameters 
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Figure 2: (a) Distribution of permutation T values and nominal T value for a simulated sample with 2% CP 
violation in the amplitude of the p + (770) resonance. Permutation T s are fitted with a GEV function and the 
nominal T value is shown as a vertical line, (b) 7] value distributions , and (c) local asymmetry significances 
are plotted on Dalitz plots for the same simulated sample. 


within their uncertainties, taking into account their correlations, and by extracting a p-value for 
each of these samplings. 

The contribution of each single event from one flavour, 7}, and from the opposite flavour %, to 
the total T value is defined as 


2 n (n 
T ‘ = 2n(n 

A permutation method is also used here to define the level of significance of each event. The 
visualisation of regions with significant asymmetries is obtained by plotting the asymmetry signif¬ 
icance in terms of each event’s contribution on Dalitz plot (Fig. 2(b)(c)). 

The practical limitation of this method is that the number of mathematical operations scales 
quadratically with the sample size. Furthermore, a significant number of permutations is required 
to get a sufficient precision of the p-value. In this analysis, graphics processing units (GPUs) [7] 
are used in the energy test calculation. The parallel architecture of GPUs reduces calculation time 
by about 95% in this analysis. 

4. Sensitivity studies 

A study of sensitivity with different types and amounts of asymmetries is required in the interpreta¬ 
tion of results. This sensitivity study is done with simplified Monte Carlo data generated according 
to the model described in Ref. [2] with the generator package Laura+-i- [8]. 

The selection efficiency is also taken into account in the sensitivity study as it varies strongly 
across phase space. The selection efficiency is obtained using full LHCb detector simulation. The 
sensitivity studies are performed using a cocktail of simulated signal and background events. The 
signal sample is generated using Laura+-i- taking into account the Dalitz plot efficiency model, 
while the background sample is generated according to the sidebands distribution. 
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Table 1: Overview of sensitivities to various CP violation scenarios. A A and A (j) denote, respectively, change 
in amplitude and phase of the resonance R. 


R (AA, A (j>) 

p —value 

(fit) 

P° 

(4%, 0°) 

3 3 +L1 

X 

10 -4 

P° 

(0%, 3°) 

1 5+ 1 - 7 

X 

10 -3 

P + 

(2%, 0°) 

5 0+ 8 - 8 

J ’ U -3.8 

X 

icr 6 

P + 

(0%, 1°) 

6 3+ 55 
°- 3 —3.3 

X 

icr 4 

P~ 

(2%, 0°) 

2 0 +L3 
z - u —0.9 

X 

icr 3 

P~ 

(0%, 1.5°) 

8 9 +22 

0>y -6.7 

X 

icr 7 


Table 2: Results for various metric parameter values. The p- values are obtained with the counting method. 
All p-values are found to be consistent with no CP violation hypothesis. 


a [GeV 2 /c 4 ] 

p —value 

0.2 

(4.6 ±0.6) x 10 -2 

0.3 

(2.6 ±0.5) x 10“ 2 

0.4 

(1.7±0.4) x 10“ 2 

0.5 

(2.1 ±0.5) x 10 -2 


A study of the metric parameter a is also performed in the sensitivity studies, a plays the role 
of bin size in a binned method. It should be larger than the resolution as otherwise asymmetries 
would be washed out by fluctuations within the scale of the metric. At the same time it should be 
small enough to resolve the regions where asymmetries lead to different densities. Within these 
two constraints, various values of a are tested and a = 0.3 GeV 2 /c 4 is chosen as the default value. 

With 100 permutations, the studies of different CP asymmetry scenarios and their sensitivities 
are reported in Table 1 . 

5. Results and conclusions 

By counting the fraction of permutation T values above the nominal T value with 1000 permuta¬ 
tions, a p -\alue of (2.6 ±0.5) x 10 -2 for default a (0.3 GeV 2 /c 4 ) is obtained. Varying the metric 
parameter results in the ^-values listed in Table 2. All results are consistent with the no CP vio¬ 
lation hypothesis. The nominal T value for data and permutation T values for no CP asymmetry 
hypothesis, and the visualisation of significant regions on Dalitz plot are shown in Fig. 3. 

In summary, a search for CP violation in the D° -a 7r - 7T + 7T 0 decay is performed using a novel 
unbinned model-independent technique, the world’s best sensitivity of this decay from a single 
experiment is achieved. The data are found to be consistent with the no CP violation hypothesis 
with a p-v alue of (2.6 ±0.5) x 10 -2 . 
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Figure 3: (a) Permutation (no CP asymmetry) T value distribution showing the fit function and the mea¬ 
sured T value as a vertical line, (b) Visualisation of local asymmetry significances. The positive (negative) 
asymmetry significance is set for the D° candidates having positive (negative) contribution to the measured 
T value, respectively. 
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